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the  Phase  Transition  Problems  * 
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Morgantown,  WV  26506 

Abstract 

In  this  paper  we  first  summarize  the  earlier  results  on  the  slow  motion  in  the 
Korteweg  theory  of  cappilarity  in  the  one-dimensional  case  and  show  some  numerical 
results.  In  the  multidimensional  case  we  discuss  the  existence  of  local  solutions  to 
the  system  of  equations  for  compressible  fluids  of  Korteweg  type. 

1  Introduction 

In  order  to  model  the  capillarity  effect  of  materials,  Korteweg  [12]  formulated  a  constitutive 
equation  for  the  Cauchy  stress  that  includes  density  gradients.  It  turns  out  that  his  theory 
is  useful  to  discuss  phase  transition  problems. 

First,  we  discuss  the  one-dimensional  isothermal  motion.  In  this  case  the  equation  we 
discuss  is  given  by 

(1.1)  utt  =  ct(ux)x  +  vuxxt  -  Suxxxx,  0  <  x  <  1,  t  >  0. 

where  it  is  the  displacement  and  uxxt  and  uxxxx  terms  represent  the  viscosity  and  the  cap¬ 
illarity  effects,  respectively.  Typical  boundary  conditions  come  from  either  a  soft  loading 
device  or  a  hard  loading  device.  Although  the  slow  motion  occurs  in  both  cases,  in  this 
note  we  discuss  the  soft  loading  case  only  for  simplicity.  The  boundary  conditions  in  this 
case  are  given  by 

u(0,  t)  —  0,  cr(ux)  d"  Wxt  c  uxxx|x=j  —  P , 

L2  Urr(0,f)=0,  U«(M)  =  0. 

The  initial  conditions  are  given  by 

(1.3)  u(x,0)  =  /(x),  ut(x,0)  =  g(x), 

where  f,g  €  Hl{ 0,1).  The  boundary  conditions  (1.2a)  show  that  the  stress  P  is  applied 
at  x  =  1.  The  boundary  conditions  (1.2b)  are  the  natural  boundary  conditions  for  the 
corresponding  variational  problem. 

*The  author  was  supported  in  part  by  Army  Grant  DAAL  03-89-G-0088. 
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In  what  follows,  we  assume  that  a  is  given  by  Fig.  1.1.  In  this  figure  (O.a*]  and 
[0",oo)  are  called  the  a-phase  and  the  ,J-phase,  respectively.  They  correspond  to  the 
different  phases  of  the  material.  The  interval  (a’,/3*)  is  called  the  spinodal  region  and 
physically  unstable.  We  denote  by  a,  8,  and  0  the  values  of  uT  at  the  intersections  of 
y  =  P  and  y  —  cr(ux)  in  the  a-phase,  the  spinodal  region,  and  the  /3-phase,  respectively. 
The  value  of  P  for  which  areas  A  and  B  are  equal  is  called  the  Maxwell  line.  We  denote 
by  0m,  and  8m  the  values  of  a,  0,  and  8,  respectively,  for  which  we  have  the  Maxwell 
line  construction. 


y 


Figure  1.1 

The  capillarity  term  was  first  introduced  by  Korteweg  [12].  Recently,  various  effects  of 
this  term  have  been  discussed.  For  example,  Serrin  [15],  [16]  reconsidered  the  Korteweg 
theory  and  has  shown  the  existence  of  steady  profile  connecting  the  a-phase  and  the  0- 
phase.  Slemrod  [17]  and  Hagan  and  Slemrod  [9]  considered  the  existence  of  travelling  wave 
solutions.  The  static  problems  concerning  the  soft  loading  case  and  the  hard  loading  case 
have  been  discussed  in  [3]  and  [4],  respectively.  The  dynamical  aspects  of  these  loading 
cases  are  discussed  in  Hattori  and  Mischaikow  [10]  and  Andrews  and  Ball  [1], 

In  Section  2  we  summarize  the  result  in  [11]  about  a  slow  motion  of  (1.1)  resembling 
the  dynamics  of  (2.3)  discussed  in  [7],  [5],  [2],  and  [8].  In  Section  3  we  show  some  numer¬ 
ical  examples  of  slow  motions.  In  Section  4  we  dsicuss  the  existence  of  local  solutions  to 
the  system  of  equations  for  two  dimensional  isothermal  motion  of  compressible  fluids  of 
Kort°w<»g  type.  The  higher  order  terms  of  density  (or  the  deformation  gradient)  in  the 
Cauchy  stress  tensor  is  not  in  general  compatible  with  the  classical  theory  of  thermody- 
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namics.  Dunn  and  Serrin  [6]  introduced  the  concept  of  interstitial  working  and  derived 
the  Cauchy  stress  tensor  compatible  with  thethermodynamics.  First,  we  summarize  their 
results  and  derive  the  system  of  equations.  Then,  we  state  the  theorem  for  the  existence 
of  local  solutions. 


2  Slow  motions  one-dimensional  case 

In  this  section  we  summerize  the  results  in  [2],  [8],  and  [11].  Multiply  (1.1)  by  uf,  integrate 
in  x  and  t,  and  then  integrate  by  parts  using  (1.2).  After  dividing  by  e,  we  have, 

(2.1)  E[u](t)  +  -f  f  i/ult(x,s)dxds  =  £[u](0), 

£  JO  Jo 

where 

(2.2)  £(«](<)=  l\^-u2t+-(W(uI)-Prc,)  +  iul)(x,t)dx. 

Jo  t.t  £  Z 

In  (2.2)  W(ur)  is  a  primitive  of  a.  For  the  remainder  of  the  paper  we  shall  assume  that  P 
=  cr(ajvr).  This  implies  that  W(ux)  —  Pux  will  be  double- well  potentials  with  equal  depth. 
For  the  sake  of  simplicity  we  shall  also  assume  that  W(ux)  —  Pux  is  given  by  (ux  —  l)2. 
The  same  conclusions  will  hold  for  more  general  non-linearities. 

Observe  that  (2.1)  is  similar  to  that  for  the  parabolic  equation 

(2.3)  I'v t  =  t2vxx  -  (v3  -  v), 
with  either  the  homogeneous  Neumann  boundary  condition 

ux(0,f)  =  0,  vx(l,t)  =  0 

or  a  Dirichlet  condition 

v(0,t)  =  a,  v(l,t)  =  b,  a,b=±  1. 


In  particular  the  energy  relation  for  (2.3)  is  given  by 


(2.4) 

1  rt  yl 

Ep[v](t)  +  -  /  J/Uf2(x,s)  dxds  =  £p[t;](0), 

£  Jo  Jo 

where 

Ep[v}{t )  =  jf  {~{v7  -  l)2  +  ^vl}{x,t)dx. 

(2.5) 

Now  we  summarize  the  results  of  the  slow  motions  for  the  parabolic  equation.  We 
assume  for  the  initial  data  of  (2.3)  that 

(2.6)  ru(x)  =  limv'(x,  0) 
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exists  as  a  limit  of  L 1  norm,  where  w  is  a  piecewise  constant  function  taking  only  the  values 
±1,  with  exactly  A  discontinuities  at  { j: t ,  ■  •  •  ,i,vj  and  we  also  assume  that  the  initial  data 
satisfy 

(2.7)  £P[t/](0)  <  Nc0  +  K2  exp(-Kje). 

Then,  we  have 

Lemma  2.1  Suppose  the  initial  data  for  (2.3)  satisfy  (2.6)  and  (2.7).  Then,  for  any  T 
satisfying  0  <  T  <  Fve*  exp(-K/e),  we  have 

(2.8)  sup  I'  \v\x,t)-vc(x,Q)\dx  <  (FG)l/2e^s+1]. 

0  <t<TJO 


Next,  we  summarize  the  results  concerning  the  slow  motions  of  (1.1).  As  the  form  of  the 
energy  relation  (2.1)  resembles  (2.4),  we  can  expect  to  draw  the  same  kind  of  conclusions 
for  (1.1).  For  this  purpose  we  rewrite  the  energy  Ec[u]  as 

Ec[u]  =  E,[u]  +  Ep[ux],  £a[u]  =  /  ^-u\(x,t)dx. 

J  0  ^6 

We  assume  that  the  initial  data  for  (1.1)  satisfy 

(2.9)  uxixi  0)  =  u£(x,  0) 
and 

(2.10)  £s[u'](0)  <  C exp(— Jf/e). 

The  condition  (2.9)  is  imposed  for  the  sake  of  simplicity.  As  long  as  ux(x,0)  satisfies  (2.6) 
and  (2.7),  with  v  being  replaced  by  ur,  the  same  conclusion  should  be  obtained. 

Lemma  2.2  Suppose  (2.9)  and  (2.10)  are  satisfied.  Then,  for  any  T  satisfying  0  <  T  < 
Fctse3  exp(K/e),  the  solution  to  (1.1)  satisfies 

(2.11)  sup  /  |<(x,  t)  -  u' (*,  0)|  dx  <  (FcGcy'2e^+1\ 

o<t<rJo 


Using  Nirenberg’s  inequality  and  Lemmas  2.1  and  2.2,  we  can  show 

Theorem  2.3  If  s  >  1,  then  the  difference  in  the  L°°  norm  between  ucx  and  vc  is  (9(e^5-1^) 
for  at  least  0  <  t  <  Fmvt3  exp(/i/e),  where  Fm  ~  min {F,FC}. 


134 


3  Numerical  examples 

We  give  a  numerical  exmaple  of  the  soft  loading  case  to  confirm  the  results  in  the  previous 
section.  We  introduce  the  transform 

p  =  J  ut(x,t)dx,  q  =  ux 
similar  to  Pego’s  [14].  Then,  (1.1)  becomes 

,  Pt  =  vPxX  ~  mxx  +  <?(q)  ~  P, 

(o.l) 

Pxx • 

The  boundary  conditions  for  p  and  q  become 

Pr(0,  t)  =  0,  p(0,f)  =  0, 

9x(0,f)  =  0»  9x(l,i)  =  0. 

For  the  initial  condition,  we  consider  the  case  when 

(3.3)  p(x,  0)  =  0,  q(x,0)  =  C/(x), 

where  C  is  a  parameter  representing  the  magnitude  of  initial  data. 

As  an  example,  we  consider  the  case  when  e  =  0.01,  v  —  1.0,  and  the  initial  data  for 
the  parabolic  equation  and  for  (3.1)  are  given,  respectively,  by 

v(z,  0)  =  C( cos  2nx  +  cos  97 rx), 

p(z,0)  =  0,  q{x,0)  =  u(z,0). 

For  C  we  gave  the  following  values: 

C  =  1.0,0.5,0.1,0.01,0.001. 


(3.2) 


One  of  the  reasons  why  we  change  the  magnitude  of  the  initial  data  is  to  see  how  this 
influences  the  metastable  states.  We  should  note  that  for  either  choice  of  C  above,  the 
conditions  (2.9)  and  (2.10)  are  not  satisfied.  Nevertheless,  when  C  =  l.U,  0.5, 0.1,  v  and  q 
have  reached  the  same  metastable  state  in  each  case.  Here,  we  show  the  numerical  results 
of  C  ~  0.1,0.01  only.  In  Figures  3.1  and  3.2  we  show  how  v  and  q  evolve  for  0  <  t  <  10 
and  then  in  Figures  3.3  and  3.4  we  show  the  profiles  of  v  and  q  at  t  =  1000.  We  use  the 
solid  lines  for  v  and  the  gray  lines  for  q.  When  these  lines  overlap,  we  see  only  the  gray 
lines.  When  C  =  0.1,  they  agree  at  least  to  10-7  at  t  =  1000  and  this  agreement  continues 
at  least  until  t  =  10000.  In  these  figures,  the  values  of  x  should  be  multiplied  by  0.01. 
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Figure  3.3.  v  and  q  at  t  =  1000  for  C  =  0.1. 


Figure  3.4.  v  and  q  at  t  =  1000  for  C  =  0.01. 
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4  Local  existence  in  multidimensional  case 


Dunn  and  Serrin  [6j  modified  the  Korteweg  theory  and  derived  the  following  set  of  equations 
for  the  conservation  of  mass,  the  balance  of  linear  momentum,  the  balance  of  energy,  and 
the  Clausius-Duhem  inequality: 


(4.1) 


pt  +  div(pu)  =  0. 

Du  A'  T 
p~=—  =  divl, 

^  Dt 

De  _  T  .. 

p—  =  T  L  —  divq  +  divw, 

Drj  q  •  (grad#) 

^  OT  +  d.vr,  +  3-il -  >  0, 


Df 

where  -  =  ft  +  u  •  V/  and 

1.  p  =  p(x,  t)  is  the  density  of  the  fluid  at  the  point  x  at  time  t, 

2.  u  =  u(x,  t)  is  the  velocity  of  flu  d, 

3.  9  =  0(x,  t)(>  0)  is  the  absolute  temperature, 

4.  e  =  e(x,  t)  is  the  specific  internal  energy  per  unit  mass, 

5.  t)  =  t/(x,  <)  is  the  specific  entropy  per  unit  mass, 

6.  T  =  T(x,  t)  is  the  Cauchy  stress  tensor, 

7.  q  =  q(x,  t)  is  the  heat  flux  vector, 

8.  L  =  gradu. 


The  main  difference  with  the  classical  thermodynamics  is  the  divw  term  and  w  is  called  the 
interstitial  work  flux  representing  spacial  interactions  of  longer  range.  They  have  proved 
that  for  a  given  Helmholtz  free  energy  ib{p,0,d),  the  following  forms  of  w  and  T 


w  =  ppipd  +  w, 

T  =  +  pd  ■  V»d  +  P2V  •  0d)I  -  pd  0  V»d 

are  compatible  with  (4.  Id).  Here,  p2x^p(p,  0,0)  is  the  pressure  and  w  is  the  “static”  portion 
of  the  interstitial  work  flux  w.  They  have  shown  that  if  the  material  possesses  a  center  of 
symmetry,  w  =  0.  In  what  follows,  we  consider  the  materials  which  posess  the  center  of 
symmetry.  They  also  have  ob;  ved  that  the  classical  forms  of  viscosity  and  conductivity 
tensors  are  compatible. 
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In  this  note  we  state  a  result  concerning  the  existence  of  a  unique  local  smooth  solution 
in  the  two-dimensional  isothermal  motion  of  the  Korteweg  type  materials  where  the  viscous 
effect  is  also  included.  The  3-dimensional  case  can  be  discussed  similarly.  In  what  follows, 
we  state  the  assumptions  on  the  Helmholtz  free  energy  and  derive  the  system  that  we  shall 
discuss.  We  assume  that  the  Helmholtz  free  energy  is  given  by 


(4.3) 


i'  =  F(p)  +  —(pi  +  p\)t 

Ip 


where  F  is  a  smooth  function  of  p  and  1/  is  a  positive  constant.  This  choice  is  to  make  the 
terms  appearing  in  (4.4)  as  simple  as  possible,  yet  to  reflect  the  effect  of  the  higher  order 
terms  of  p.  1 

With  the  choice  the  Helmholtz  free  energy  given  in  (4.3)  and  with  A  =  —  -p,  the  system 

o 


then  becomes 
(4.4) 


Pt  +  (pu)x  +  (/”>)»  =  0, 

(pu) t  +  {pU-2)x  +  {Pnv)y  ~  (TU)x  +  (^12 )y i 

(pv) t  +  ( PUV)X  +  (pv2)y  =  (T2  l)x  +  (T2  2/y, 
where  u  and  v  are  the  x  and  y  component  of  velocity  and 

Tu 
T2  2 


(4.5) 

T  —  (T * 

\t2X 

=  {-P  + 

(4.6) 

and 

V  = 

(4.7) 

pI 

PxPy 


PxPy 


+  v, 


P  =  P2F\p), 


Vu  Vl2 
V2i  v22 


=  p{(gradu)  +  (gradu)T  -  -(divu)I}. 


Here,  I  is  the  unit  rank-two  tensor  and  a  superscript  T  denotes  the  transpose  of  a  tensor. 
Since  we  discuss  the  existence  of  a  local  solution,  we  do  not  need  the  monotonicity  of  the 
pressure  on  p.  Further  computation  simplifies  the  divT  term 


(4.8)  divT  =  — Vp  +  i/pV(Ap)  -f  divV. 

We  discuss  the  local  existence  for  the  pure  initial  value  problem  of  (4.4)  with  the  initial 
data  given  by 

(4.9)  (p,tt,v)(x,2/,0)  =  {p0,u0,v0){x,y). 

Another  reasonable  choice  ic  to  change  the  last  term  in  (4.3)  with  ^(pl  +  Py)-  Although  this  choice 
may  be  physically  more  realistic,  mathematically  it  is  more  cumbersome  to  handle.  For  example,  the 
expression  for  divT  is  very  complicated.  Therefore,  we  do  not  discuss  this  case  (See  (4.8)). 
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We  assume  that  the  initial  data  satisfy 

(4.10)  (po  —  Po,  ^0,  t’o)  €  Hk(R2),  p0  >  S  >  0, 

where  k  >  4  and  po  >  0  is  a  positive  constant.  Denote  by  |)  •  ||  =  ||  ■  j]0  the  I2  norm  and 
by  |j  •  Ijk  the  k- th  order  Sobolev  norm.  Set 

M2.t=  sup  (u«KOIIJ  +  l|Vp(0ll2)  +  [T  (l|Vu(<)||!  +  ||Vu(i)f )  it 

0 <t<T  V  '  Jo  v  ' 

and 

MI  =  £ 

bi<* 

where  w  =  (p,u,v).  The  main  result  is  stated  as  follows. 

Theorem  4.1  For  any  initial  data  (po,uQ,vo)  such  that  p0  >  S  >  0  and  ( p0  —  po,u0,i'o)  G 
Hk(R 2)  ( k  >  4)  where  po  >  0  is  a  constant,  there  exists  a  T  >  0  such  that  in  t  6 
[0,T],  the  Cauchy  problem  (4.4),  (4.9)  has  a  unique  solution  ( p,u,v )  such  that  p  —  p0  € 
L°°([ 0,  T];  Hk+l(R2))  and  (u,  v)  G  L°°([0,  T\;  Hk(R2))  and 

Mtl<ck\\w0\\l  +  \\P0\\l+l. 

Since  the  linearized  problem  of  (4.4),  (4.9)  is  not  of  any  classical  type,  the  existence  of 
solutions  is  not  known  even  for  the  linearzed  problem.  We  prove  the  existence  of  solutions 
for  the  linearized  problem  by  establishing  an  energy  estimate  for  the  dual  problem  and 
then  using  the  dual  argument. 
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